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Abstract 

The flow of Ree-Eyring and Casson non-Newtonian fluids is investigated 
using a variational principle to optimize the total stress. The variationally- 
obtained solutions are compared to the analytical solutions derived from the 
Weissenberg-Rabinowitsch-Mooney equation and the results are found to be 
identical within acceptable numerical errors and modeling approximations. 
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1 Introduction 

Many analytical and numerical methods have been invented and used, especially in 
the recent decades, for modeling, simulating and analyzing the flow of various types 
of fluid under diverse flow conditions in conduits of different shapes and geometries 
and with different mechanical properties. These methods include the application of 
first principles of mechanics, the use of Navier-Stokes equation, the application of 
lubrication approximation, the use of Weissenberg-Rabinowitsch- Mooney relation, 
stochastic methods, and finite element and spectral methods [1-9]. 

Recently, we proposed [10] the use of a new method which is based on applying 
the Euler-Lagrange variational principle to optimize the total stress of fluid in 
the flow conduit for the purpose of obtaining relations for the flow of generalized 
Newtonian fluids. Although the method is general, it lends itself more easily to the 
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flow in rigid straight uniform conduits with circularly-shaped cross sections as the 
flow in these cases is essentially unidirectional. 

There, the method was demonstrated by using a straight cylindrical tube ge- 
ometry with six fluid rheological models: Newtonian, power law, Bingham plastic, 
Herschel-Bulkley, Carreau and Cross. The proposed method can be applied ana- 
lytically, when the variational equation can be solved in all its stages by analytical 
means, as well as numerically in part when some stages are difficult or impossible 
to solve analytically and hence standard numerical methods, like numerical inte- 
gration by Simpson’s rule and solving implicit equations by bisection methods, can 
be used. It was also demonstrated that although the variational principle in its 
proposed form is restricted to fluids, and hence it does not include viscoplastic 
materials which behave like solids prior to their yield point, the method is a good 
approximation when applied to the flow of viscoplastic materials of low yield stress 
value since the effect of the non- yielded plug region on the total flow is minimal. 

In a later study [11] we provided more support to the variational method by 
successful application to the flow in a plane long thin slit geometry using two 
rheological models: Newtonian and power law where the well-known analytical 
solutions for these fluids were obtained from the application of the variational 
principle. 

In the present study, we further investigate and validate the variational method 
through its application to the flow of two non-Newtonian fluid models, Ree-Eyring 
and Casson, in rigid straight uniform pipes with circular cross sectional shape. The 
first of these is a pseudoplastic model which is normally used to describe the flow of 
shear thinning fluids, while the second is a viscoplastic model used to describe the 
rheology of yield stress materials. The method is applied analytically in part and 
numerically in other parts. The investigation related to the Casson model provides 
further support to the fact that the variational method is a good approximation 
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for the viscoplastic materials characterized by a low value of yield tress. 

2 Variational Method 

Before outlining the variational method and its application to the flow of Ree- 
Eyring and Casson fluids, we state the assumptions and approximations that 
have been considered in the model development. In this investigation, we as- 
sume a laminar, incompressible, steady, pressure- driven, fully-developed flow of 
time-independent, purely- viscous fluids that can be described by the generalized 
Newtonian fluid model which is given by 


t = H 7 ( 1 ) 

where r is the shear stress, // is the generalized Newtonian viscosity, and 7 is the 
rate of shear strain. As for the conduit, we assume a rigid straight tube with a 
uniform cross section (i.e. constant shape and area along the axial dimension) of 
circular shape. 

The variational method for obtaining the flow relations for generalized New- 
tonian fluids is based on the minimization of the total stress in the flow conduit 
using the Euler-Lagrange variational principle to obtain the rate of shear strain. 
This is then followed by obtaining the flow velocity profile by integrating the rate 
of strain with respect to the tube radius, which is followed by integrating the flow 
velocity profile with respect to the cross sectional area to obtain the volumetric flow 
rate. Although the first step in this method is based on analytical manipulation of 
the derived variational principle, the subsequent stages can be manipulated either 
analytically or numerically. This allows more flexibility in the application of the 
variational method as will be vividly demonstrated in the forthcoming paragraphs 
and sections. 
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From the basic principles and definitions of the fluid mechanics of generalized 
Newtonian fluids, the total stress arising from the flow profile of a fluid flowing in 
a tube of radius R is given by 


n = 




07) dr 



da 
dr 



dr 


( 2 ) 


where r t is the total shear stress, r m and t w are the shear stress at the tube center 
and wall respectively, and r is the radius at an arbitrary point on the tube cross 
section. 

The total stress, as given by Equation 2, can be optimized by applying the 
Euler-Lagrange variational principle which, in one of its diverse forms, is given by 
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where the symbols in the last equation correspond to the symbols of the stated 
flow problem as follow 
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Considering the fact that for the considered flow systems 


dy d'y 


G 


( 6 ) 


where G is a constant, it can be shown that Equation 5 can be reduced to two 


independent variational forms 
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In the following two subsections we use the second form where we outline the ap- 
plication of the variational method, as summarized in Equation 8, to the Ree-Eyring 
and Casson fluids. We then validate the variationally-obtained flow solutions 
by comparing them to the analytically-derived solutions from the Weissenberg- 
Rabinowitsch-Mooney (WRM) equation. 

2.1 Ree-Eyring Fluids 

For the Ree-Eyring fluids, the rheological constitutive relation that correlates the 
shear stress to the rate of shear strain is given by 


where /r 0 is the viscosity at vanishing rate of shear strain, and r c is a characteristic 
shear stress. Hence, the generalized Newtonian viscosity is given by 



(9) 
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On substituting // from the last relation into Equation 8 we obtain 



( 11 ) 


On integrating once we get 
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t c asinh 
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( 12 ) 
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dr 


where A is a constant. On separating the two variables and integrating again we 
obtain 


where the constant of integration C is absorbed within the integral on the left hand 
side. The left hand side of Equation 13 is not a standard integral that can be easily, 
if possible at all, to integrate analytically using standard integration methods or 
can be found in tables of integrals. On using a computer algebra system, the 
following expression was obtained 


where Li 2 is the polylogarithm function. This expression, when evaluated, produces 
significant errors as it diverges with decreasing r due apparently to accumulated 
numerical errors from the logarithmic and polylogarithmic functions. 

To solve this problem we used a numerical integration procedure to evaluate 
this integral, and hence obtain A, numerically using the boundary condition at the 
tube wall, that is 



(13) 
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where t w is the shear stress at the tube wall which is defined as the ratio of the force 
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normal to the tube cross section, F±, to the area of the luminal surface parallel to 
this force, A», that is 

F± nR 2 Ap RAp 

T - = t, = m = (16) 

where R and L are the tube radius and length respectively, and A p is the pressure 
drop across the tube. 

This was then followed by obtaining 7 as a function of r using a bisection 
numerical solver in conjunction with a numerical integration procedure based on 
Equation 13. Due to the fact that the constant of integration C is absorbed in 
the left hand side and a numerical integration procedure was used rather than an 
analytical evaluation of the integral on the left hand side of Equation 13, there was 
no need for an analytical evaluation of this constant using the boundary condition 
at the tube center, i.e. 

7 (r = 0) = 0 (17) 

The numerically-obtained 7 was then integrated numerically with respect to r to 
obtain the flow velocity as a function of r. The flow velocity profile was then inte- 
grated numerically with respect to the cross sectional area to obtain the volumetric 
flow rate, as outlined previously. 


2.2 Casson Fluids 


For the Casson fluids, the rheological constitutive relation is given by 


t = 


(Vr ) 1/2 + T o ' 2 


(18) 


where k is a viscosity consistency coefficient, and ro is the yield stress. Hence 


7 


r 



(fry) 


1/2 , _l/2 
■+■ r 0 


7 


( 19 ) 


On substituting // from the last equation into the main variational relation, as given 
by Equation 8, we obtain 
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On integrating once we get 
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where A is a constant, that is 
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On separating the two variables and integrating again we obtain 
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(23) 


k'y + 4 (fcT 0 7) 1,/2 + r 0 In (7) = Ar + C (24) 

where C is another constant. 

Now, we have two boundary conditions one at the tube center and the other at 
the tube wall. As for the first we have 


7 (r = 0) = 0 


(25) 


8 


As the solution should be finite for all values of r including r = 0, we should have 


C = 0 


(26) 


As for the second boundary condition we have 


k^ w + 4 (kT 0 y w ) 1/2 + r n 111 (7 W ) = AR 


(27) 


where the rate of shear strain at the tube wall, 7^, is given by 



1/2 


k 


(28) 


with the shear stress at the tube wall, t w , being given by Equation 16. ffence 


Equation 24, which defines 7 implicitly in terms of r, is then used in conjunction 
with a numerical bisection method to find 7 as a function of r. The numerically 
found 7 is then integrated numerically with respect to r to obtain the flow velocity 
profile which is then integrated numerically with respect to the cross sectional area 
to find the volumetric flow rate. 


The method, as described in the last section, was implemented and applied to 
the Ree-Eyring and Casson fluids over an extended range of the tube and fluid 
parameters, and the results were examined and analyzed. In all cases we obtained 
very good agreements between the variational solutions and the WRM analytical 
solutions; a sample of the investigated cases are presented in Figures 1 and 2. 

The difference between the two solutions for the Ree-Eyring fluids can be largely 


A _ kj w + 4 ( fcr 0 7^) 1/2 + t 0 In (j w ) 


(29) 


3 Results and Analysis 
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explained by numerical errors arising from the extensive use of numerical integra- 
tion and bisection solvers in several stages of the method implementation. These 
errors propagate and amplify on passing from one stage to the next. However, 
almost in all cases that have been investigated, the maximum and the average of 
the percentage relative difference between the variational and the WRM analytical 
solutions do not exceed 1%. 

As for the Casson fluids, as well as the above mentioned accumulated and mag- 
nified errors from the persistent use of numerical integration and numerical solvers 
in multiple sequential stages, there is a more fundamental reason for the deviation 
between the two solutions, that is the variational method is strictly applicable to 
real fluids with no yield stress and hence the use of this method for the viscoplastic 
materials is just an approximation which is usually a good one when the value of 
the yield stress is low. Hence, in Figure 2 we can see an obvious and logical trend 
that is as the yield stress value increases, the deviation between the variational and 
the WRM analytical solutions increases as it should be because the approximation 
in using the variational method is worsened as the value of the yield stress rises. 
As the yield stress increases, the non-yiclded plug region in the middle of the tube, 
which is not strictly subject to the variational formulation, increases in size and 
hence the gap between the two solutions is widened. 
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(a) (b) 



(c) (d) 

Figure 1: Comparing the WRM analytical solution, as given by Equation 35, to the 
variational solution for the flow of Ree-Eyring fluids in circular straight uniform 
rigid tubes with (a) R = 0.02 m, L = 0.5 m, /i 0 = 0.01 Pa.s, and r c = 500 Pa, 
(b) R = 0.01 m, L = 0.7 m, Ho = 0.005 Pa.s, and r c = 750 Pa, (c) R = 0.003 m, 
L = 0.05 m, /i 0 = 0.02 Pa.s, and r c = 300 Pa, and (d) R = 0.006 m, L = 0.08 m, 
Ho = 0.015 Pa.s, and r c = 400 Pa. In all four sub-figures, the vertical axis represents 
the volumetric flow rate, Q, in m 3 .s -1 while the horizontal axis represents the 
pressure drop, A p, in Pa. The average percentage relative difference between the 
variational and the WRM analytical solutions for these cases are about 0.67%, 
0.65%, 0.77%, and 0.75% respectively. 
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(a) (b) 



(c) (d) 


Figure 2: Comparing the WRM analytical solution, as given by Equation 39, to 
the variational solution for the flow of Casson fluids in circular straight uniform 
rigid tubes with (a) R = 0.01 m, L = 0.1 m, k = 0.005 Pa.s, and r Q = 0.1 Pa, 
(b) R = 0.05 m, L = 0.5 m, k = 0.07 Pa.s, and r 0 = 0.25 Pa, (c) R = 0.05 m, 
L = 0.5 m, k = 0.07 Pa.s, and r D = 0.75 Pa, and (d) R = 0.01 m, L = 0.1 m, 
k = 0.005 Pa.s, and r D = 1.0 Pa. In all four sub-figures, the vertical axis represents 
the volumetric flow rate, Q , in m 3 .s -1 while the horizontal axis represents the 
pressure drop, A p, in Pa. The average percentage relative difference between the 
variational and the WRM analytical solutions for these cases are 0.55%, 1.03%, 
2.26%, and 4.07% respectively. 
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4 Conclusions 


In this paper, we demonstrated the use of the variational method that optimizes 
the total stress of the fluid in the flow conduit to obtain flow relations corre- 
lating the volumetric flow rate to the pressure drop for the flow of two non- 
Newtonian fluid models, namely Ree-Eyring and Casson, in rigid straight uniform 
tubes with circular cross sectional shape. The variational method was closely exam- 
ined and analyzed using extensive ranges of fluid and conduit parameters and the 
results were compared to the analytically-obtained solutions from the Weissenberg- 
Rabinowitsch-Mooney equation. In all cases, good agreement was observed between 
the variational and the WRM analytical solutions within anticipated and justifiable 
numerical errors and modeling approximations. 

Although the variational principle in its current formulation strictly applies 
only to non-viscoplastic fluids, the analysis shows that the variational method 
is a good approximation for the viscoplastic fluids with low value of yield stress 
as demonstrated by several examples from the Casson fluids. An obvious logical 
trend was observed where the gap between the variational and the WRM analytical 
solutions was enlarged as the value of the yield stress of the viscoplastic fluids was 
increased. 

Overall, this investigation adds more support to the previous investigations 
[10, 11] related to the application of the Euler-Lagrange variational principle to 
the flow of generalized Newtonian fluids through rigid straight uniform pipes with 
circular cross sections and slits of a plane long thin geometry. Further investigations 
are planned for the future to extend the application of this method to other types 
of fluid, flow condition and flow conduit to verify the validity of the variational 
principle and the applicability of the derived method to the flow dynamic systems 
in general. As part of this investigation, the volumetric flow rate relations for the 
Ree-Eyring and Casson fluids in circular tubes were derived using the well-known 
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Weissenberg-Rabinowitsch- Mooney equation. 

Nomenclature 

7 rate of shear strain 

r y w rate of shear strain at tube wall 

H viscosity of generalized Newtonian fluids 

/io viscosity at zero shear rate in Ree-Eyring model 

r shear stress 

To yield stress in Casson model 

r c characteristic shear stress in Ree-Eyring model 

r m shear stress at tube center 

Tt total shear stress 

t w shear stress at tube wall 

A\\ area of tube luminal surface 

F± force normal to tube cross section 

k viscosity consistency coefficient in Casson model 

L tube length 

Li 2 polylogarithm function 

p pressure 

A p pressure drop 

Q volumetric flow rate 

r radius 

R tube radius 
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A Deriving Flow Relations for Ree-Eyring and 
Casson Fluids 

In this appendix we derive analytical expressions correlating the volumetric flow 
rate to the pressure drop for the flow of Ree-Eyring and Casson fluids in rigid 
straight uniform tubes with circular cross sections which we could not find in the 
literature. For this purpose we use the Weissenberg-Rabinowitsch-Mooney relation 
which is given by [12] 


Q = nI j- [ T 2 ^dr (30) 

Tw JO 

where Q is the volumetric flow rate, r is the shear stress, 7 is the rate of shear 
strain, and r w is the shear stress at the tube wall which is given by 


7~W 


RAp 

2 L 


(31) 


where R and L are the tube radius and length respectively, and A p is the pressure 
drop across the tube. 


For the Ree-Eyring fluids we have 

r = t c asinh (^~^j (32) 

where po is the viscosity at vanishing rate of shear strain, and r c is a characteristic 
shear stress. Therefore 
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Inserting this into the Weissenberg-Rabinowitsch-Mooney relation we obtain 
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